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Abstract 

We have constructed a Doppler-shifter-type pulsed ultra-cold neutron (UCN) source at the Ma¬ 
terials and Life Science Experiment Facility (MLF) of the Japan Proton Accelerator Research 
Complex (J-PARC). Very-cold neutrons (VCNs) with 136-m/s velocity in a neutron beam supplied 
by a pulsed neutron source are decelerated by reflection on a m=10 wide-band multilayer mirror, 
yielding pulsed UCN. The mirror is fixed to the tip of a 2,000-rpm rotating arm moving with 
68-m/s velocity in the same direction as the VCN. The repetition frequency of the pulsed UCN is 
8.33 Hz and the time width of the pulse at production is 4.4 ms. In order to increase the UCN flux, 
a supermirror guide, wide-band monochromatic mirrors, focus guides, and a UCN extraction guide 
have been newly installed or improved. The 1 MW-equivalent count rate of the output neutrons 
with longitudinal wavelengths longer than 58 nm is 1.6 x 10^ cps, while that of the true UCNs 
is 80 cps. The spatial density at production is 1.4 UCN/cm^. This new UCN source enables us 
to research and develop apparatuses necessary for the investigation of the neutron electric dipole 
moment (nEDM). 


1 


1. INTRODUCTION 


Neutrons with kinetic energies lower than the Fermi pseudo-potential of Ni are called 
“ultra-cold neutrons (UCNs).” The typical value of this potential is 243 neV, while the de 
Broglie wavelength, A, of a UCN is more than 58 nm. A UCN experiences a potential (which 
represents the neutron-nucleus interaction spatially averaged over millions of nuclei) having 
a form of a step function which takes a nonzero value within the space where the UCN 
wave function overlaps the potential. The incident wave-packets colliding with the averaged 
potential barrier that is higher than the kinetic energy of the neutron penetrate into the 
classically forbidden region and hence some UCNs are scattered by nuclei inside the barrier, 
and then the wave-packets are almost fully reflected. The magnetic dipole moment of a 
neutron in a magnetic field has a potential energy that is expressed by the scalar product 
of the neutron magnetic moment and the magnetic flux density B, that is, —/x • B. 
The absolute value corresponds to 60 neV at a magnetic flux density of 1.0 T. It means 
that polarized UCNs with velocities less than 3.4 m/s cannot overcome the magnetic held. 
As results of these features, UCNs can be stored in a magnetic bottle or a vessel with an 
interior surface which has a high material potential, and therefore are useful for long-term 




storage experiments snch as searches for the nentron electric dipole moment (nEDM) 

^ or measurement of neutron lifetime { 4 , 5]. The non-zero nEDM implies the breaking of 
time-reversal invariance and its magnitude limits the predictions of some new physics models 
beyond the standard model of elementary particles. The neutron lifetime determines the 
mass fraction of ^He in the early Universe and the improvement of the accuracy of the 


lifetime verihes the consis 


and satellite surveys js, 7|. In addition, a UCN is useful probe in studies of a region far 
from the material surface, such as the behavior of the potential or the phonon dispersion 


ency between the fractions resulted from lifetime measurements 


relations 


8 |. 


The most well-known UCN source is Steyerl neutron turbine at Institut Laue-Langevin 
(ILL) j^. Very cold nentrons (VCNs) generated in the reactor are guided by a slightly 
curved vertical neutron guide and decelerated below 50 m/s by the gravitational held, while 
thermal neutrons are hltered out by the curved guide. The VCNs are guided into a vacunm 
chamber in which a turbine with 690 cylindrically curved blades rotates with a blade velocity 
of 25 m/s. The blade surface is coated with Ni and the neutrons are multiply rehected by 
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the moving blades. Finally, the VCNs are decelerated below 6.2 m/s by the Doppler e 


A neutron Doppler shifter is a similar apparatus to the abovementioned turbine 
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ec 


12 |. 


In this apparatus, a small number of neutron mirrors move backwards at half the velocity 
of the target neutrons, and the neutrons are decelerated to the UCN region by a single 
reflection only. The Doppler shifter is an effective device for decelerating pulsed neutrons 
with large velocity components in the guide direction. The simple mechanism allows us 
to estimate the UCN flux precisely by a particle tracking simulation. In the simplest case 
of the reflection from a rotating small flat reflector, there is a strong correlation between 
the velocities and recoil angles of decelerated neutrons. The neutrons have a short-pulse 
time structure and therefore are useful for a time-of-flight (TOF) spectrometer 13|. By 
using a shutter beside the source point whose open and close operation is synchronized with 
UCN production, the Doppler shifter can also be used for UCN storage experiments as if by 
using a steady source at the peak density of UCN [lO, To the best of our knowledge, 
the first pulsed UCN source produced by a Doppler shifter was reported by Dombeck et 


al. 


HIn, 


13|; this Doppler shifter was equipped with a package of Thermica crystals on 


the tip of a rotating arm and decelerated 400-m/s neutrons. 

Another Doppler shifter has since been developed for research and development (R&D) at 
the BL05 NOP beamline of the Materials and Life Science Experiment Facility (MLF) at the 
Japan Proton Accelerator Research Complex (J-PARC) jlJj.Il5l|. Neutron supermirrors 


16 


18| were selected for the reflector of our Doppler shifter because the reflectivity wavelength 


band is controllable freely. The peak phase space density of 136-m/s VCN beam at J-PARC 
(in Ref. 19|) was 31 times greater than that of the UCNs at ILL PF2 (in Ref. {q]), hence the 
production of extremely high density UCN was expected. The demonstration experiment of 
UCN production was carried out at a beam power of 120 kW. As a result, a UCN count of 
1.8 cps was expected for a forthcoming 1-MW beam operation 1^. The low count rate was 
due to a strong collimation of the injected VCNs and a small solid angle of UCN extraction. 

We have proposed a nEDM search with a spallation UCN source (J-PARC P33) 20|. 
Toward the realization of the plan, the Doppler shifter was used for the R&D of an impor¬ 


tant apparatus for our nEDM search, which we named “UCN rebuncher’" 


2lj. The UCN 


rebuncher is a neutron accelerator, which accelerates or decelerates pulsed UCNs suitably 


with the combination between the gradient magnetic field 


22 | and the adiabatic fast pas¬ 


sage spin flipper |23|, rotates the UCN distribution in velocity-position phase space and then 
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time-focuses the UCNs on any position. The UCN rebuncher is necessary to improve the 
storage density of UCN at the nEDM experiment with the spallation source of low repetition 
frequency like J-PARC. Because of such a mechanism of time-focus, the pulsed UCN pro¬ 
duced by the Doppler shifter is suitable for the performance test of the UCN rebuncher. In 
the test, UCNs extracted from the Doppler shifter were transported to a detector through 
Ni-coated UCN guides passing through the UCN rebuncher. The focal length of the re¬ 
buncher was set to 3.8 m. In consideration of an arrangement of the rebuncher at BL05, 
glass guides with total length of 5.6 m were produced. As a result of a UCN transport 
experiment, the 1 MW-equivalent total count rate at the end of the UCN guide was 0.1 cps 
due to the low transmittance of the guide, while that of the environmental background was 
0.03 cps. In a time-focus experiment, unnecessary UCN pulses have to be interrupted by 
a synchronized shutter near the Doppler shifter in order to prevent frame overlap of the 
TOP and detect 3-4-m/s UCNs. In that case the average flux of UCN is greatly thinned. 
Therefore it was difficult to obtain the TOP spectrum with the shutter within a practical 
time span. In order to perform the UCN rebuncher experiment at J-PARC/MLF BL05, it 
was necessary to increase the peak flux of the UCN pulse at least 10 times. 

In order to increase the UCN flux, the incident VCNs should be transported into the 
Doppler shifter more efficiently. Therefore, we have newly installed supermirror guides, 
wide-band monochromatic mirrors, and focus guides in the neutron beamline. Further, the 
produced UCNs should also be extracted without loss. Thus, a Ni-coated UCN extraction 
guide has been inserted into the Doppler shifter. The outcomes of these upgrades have 
been evaluated at beampowers of 200-300 kW. In this paper, we describe the details of the 
construction of our upgraded UCN source and its performance. At present, the improved 
UCN rebuncher has been developing. The performance test will be carried out by using the 
improved Doppler shifter. 


2. PRINCIPLES OF UCN PRODUCTION USING DOPPLER SHIFTER 

In this section, we describe the principles of UCN production using a Doppler shifter. 
As neutrons faster than UCNs can be reflected by Bragg scattering, multilayer-structured 
materials such as mica crystals or artihcial supermirrors can be applied to the Doppler shifter 
reflector. Hereafter, we call the reflector a “Doppler mirror”. 
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A neutron is reflected if it satisfies the Bragg condition in the frame of the moving Doppler 
mirror. The velocity vector, Vr, of the reflected neutron in the laboratory frame is given by 

Vr = V + 2\{v - Vm) ■ n\n, (1) 

where v is the velocity vector of the incident neutron, Vm is the velocity vector of the mirror 
surface at the reflection point of the neutron, and n is the normal vector at the reflection 
point. In the laboratory frame, the incident neutron receives a momentum stroke in the 
direction n, as shown in equation ([T]) The wavelength of the incident neutron is shifted to 
a shorter or longer value, as in the Doppler effect of a sound wave. This is the reason why 
such devices are called Doppler shifters. 

If Vm is close to v/2, the reflected neutron is decelerated close to zero velocity. The 
decelerated neutron is regarded as a UCN if the magnitude of Vr is less than 6.8 m/s. Thus, 
a mechanism that can move a mirror in the same direction as and at the half-speed of the 
incident neutron is required, in order to form a Doppler shifter that decelerates an incident 
neutron to a UCN through a single reflection. In the apparatus described here, continuous 
motion of the mirror is realized using a rotating arm, where the mirror is attached to the 
arm tip. 

Next, we formulate the velocity of the decelerated neutron. We define the beam direction 
as along the z-axis, the rotation axis of the arm is the x-axis, and the axis perpendicular to 
the X- and z-a.xes is the y-a.xis, as shown in FIG. [T^. The incident neutrons travel in the 
positive z direction. The Doppler mirror is positioned on the arm and the mirror surface 
is parallel to the x-axis and the direction of the arm. We dehne the point where the beam 
crosses the Doppler mirror when the mirror surface comes to overlap the y-axis as the pseudo 


source point. Then, the x, and ^ components of the velocity of a 

Vy, Vz-, respectively, are given by 

reflected neutron, v^, 


(2) 

Vy ~h 0^-j 

(3) 

= 2Vm - V + 3Vm0‘^, 

(4) 


where 6^ and 6y are the respective incident angles of the neutron in the x- and ^/-planes, 
V and Vm are the magnitudes of v and Vm, respectively, and 6 is the rotation angle of the 
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mirror against the |/-axis. In these equations, we assume that Vm is close to n/2 and 5x, Sy, 
and 6 are small and ignore terms of order higher than the third order of their products. 


(a) (b) 




FIG. 1: (a) Schematic diagram of Doppler shifter deceleration; (b) illustration of velocity phase 
space distribution of decelerated neutrons. 


In this case, the velocities of the decelerated neutrons are distributed on a parabola in 
the Vy-Vz velocity space plane, as shown in FIG. dJo. Then, is expressed as 

3 . ^2 


Vz = 


AVr. 


'{Vy 2,6yVi-n') (n ‘2Vm)- 


( 5 ) 


Suppose that v = 136 m/s, Vm = 68 m/s, 6^ = 0, Sy = 0, the arm length is 325 mm and, 
accordingly, the angular velocity of the rotation is 2007r/3 rad/s. Then, the condition for 
UCN production, < 6.8 m/s, is expressed as |0| < 50 mrad, from equations ([3]) 

and (jl)). The angular width of 100 mrad can be converted to a time width of 480 fis. 

UCNs are produced in the vicinity of the pseudo source point. As an increase of 6 in Vy 
can be canceled by non-zero 6y, as shown in equation ([3]), the actual width of 9 for UCN 
production is wider for a divergent beam with a width of 6y. The estimated results yielded 
by a 3D Monte Carlo simulation based on ray tracing and incorporating the geometry of 
our Doppler shifter, which we refer to as the “UCN simulation” in this paper, indicate that 
99.7% of the extractable UCNs are produced within an angular width of 170 mrad from an 
incident VCN beam with velocity divergence within ±90 mrad. Therefore, the time width 
of the extractable UCN production is 0.9 ms. 

As the produced UCNs are extracted in the direction perpendicular to the orbit of the 
mirror, the time width of the produced UCN pulse is roughly identical to the mirror size. 
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This means that the maximum time width of the UCN pulse with 6.0-m/s velocity in the 
extraction direction is 5.0 ms, for the sample case of extraction in the upper direction with a 
30-mm high mirror. If UCNs are extracted along the |/-axis in the UCN simulation, the full 
width at half maximum (FWHM) of the UCN pulse time structure for a longitudinal velocity 
of approximately 6 m/s is 2.9 ms at production, because of the incident beam divergence. 
The FWHM of the UCN pulse with a longitudinal velocity of more than 3.0 m/s is 4.4 ms. 
Therefore, the Doppler shifter with the design described above can generate sharply pulsed 
UCN of 4.4-ms width at production. 


EXPERIMENTAL APPARATUS 


At the J-PARC MLF, pulsed proton beams are injected from the 3-GeV rapid cycling 
synchrotron ring (RCS) to the MLF source. Neutrons are then produced by spallation 
reactions and cooled by liquid hydrogen moderators. The repetition frequency of the neutron 
production is 25 Hz. The neutron production pauses for 360 ms every 3 or 6 s, because proton 
beams are supplied to the 50-GeV main ring (MR) when particle physics experiments in the 
Neutrino Experimental Facility or the Hadron Experimental Facility are being performed. 

At the BL05 port, an incident neutron beam is transported by m = 2 supermirror guides, 
where m is the critical angle normalized to that of Ni. The beam is then divided into three 
branches using supermirror bending guides: the polarized, unpolarized, and low-divergence 
beam branches, which have been discussed previously in jl4j . [iSj. The Doppler shifter is 
positioned on the unpolarized beam branch as shown in FIG. [21 where the neutrons are 
guided upwards at 2.6° by a m = 3 arc-shaped supermirror bender of 100-m radius and 4.0-m 


length 


14l |. The bending guide has a cross section of 40 mm (width) x 50 mm (height), which 


is separated into hve channel structures of 40 mm (width) x 10 mm (height). Fast neutrons 
escape from the guide at the curve and only cold neutrons are transported. Neutrons 
transported through the bender pass m = 3 supermirror guides to reach the exit of the 
concrete shielding surrounding the apparatus. These guides have 4.2-m length and the 
internal cross section is 40 mm (width) x 50 mm (height). The gap between the bender and 
guide is 235 mm wide, and is designed to be as short as possible to transport VGN with 
large divergence. Note that a B 4 G collimator was used instead of the supermirror guide tube 
when the demonstration experiment discussed in this paper was performed 19 1 . 


7 






FIG. 2: Layout of supermirror guides and Doppler shifter at J-PARC MLF BL05 with distance 
from the moderator. 


We measured the fluxes of the unpolarized branch with the supermirror guides or 
the B 4 C collimators in terms of wavelength distributions using a neutron beam monitor, 
MNH10/4.2F 2^, with a collimator of 010 mm. The detector was mounted on the unpo¬ 
larized branch exit at the 16.2-m position shown in FIG. [2l The neutron pulse interval at 
MLF is usually 40 ms, so neutrons with velocities less than 400 m/s cannot normally be 
distinguished at the unpolarized branch. Therefore, neutrons produced by the last beam 
pulse immediately before the beam supply to the MR were selected from the data to pre¬ 
vent frame overlap with slower neutrons, and the TOFs from their production were then 
converted to the corresponding wavelengths. The spectra are shown in FIG. [H which were 
calculated using detection efficiency compensation and normalized to the designed beam- 
power of 1 MW. 


3.1. Supermirror guides 


It is possible to £t the flux data sets using a Boltzmann distribution under certain con¬ 
ditions, where 


0(A) = $ 




{2TTmkBf^ """P V T 




— aA I -|- Ch 


( 6 ) 


In this equation, 0(A) is the flux distribution, <F is the total flux, h is the Planck constant, 
m is the neutron mass, ks is the Boltzmann constant, T is the neutron temperature, a is 
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FIG. 3: Neutron fluxes of unpolarized branch at J-PARC MLF BL05. The closed and open 
circles indicate the fluxes obtained with the m = 3 supermirror guides and the B 4 C collimator, 
respectively. 


the attenuation coefficient (such as absorption by the A1 windows), and Cb is the neutron 
environmental background. This expression can be applied to cases where the neutron 
divergence is independent of wavelength, such that the neutrons are perfectly reflected or 
absorbed by a neutron guide or collimator. The flux is primarily restricted by the solid angle 
of the guide entrance or the collimator exit. In the case of BL05, the m = 2 supermirror 
guides with incident angle less than the critical angle have more than 90% reflectivity; hence, 
this condition is roughly satished for low-divergence VCNs. 

In a divergence measurement with a 1-mm-diameter pinhole and a position-sensitive 


detector 


25 


26j . the FWHM of the x'- and ^/'-divergence of VCNs with wavelengths longer 


than 2.0 nm (198 m/s) were found to be 50 and 54 mrad respectively, where the x'- and 
//'-axes were the horizontal and vertical directions on the plane parallel to the unpolarized 
branch window. The divergence of the 99% VCNs was within ±60 mrad in each direction. 
Because the critical velocity of the m = 2 supermirror at zero incident angle is 13.6 m/s, 
a VCN with a glancing angle of 60 mrad and v < 13.6/ tan(0.06) = 226 m/s (1.74 nm) 
is perfectly reflected by the surface. Accordingly, the flux attenuation for such a VCN is 
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almost independent of its wavelength in the transport, excluding the case of neutron capture 
in Al. Therefore, the spectra over 2.0-nm wavelength shown in FIG. |3]are described by the 
Boltzmann distribution of equation ([2]). In the htting, we assumed the existence of a constant 
Gft = 1000 neutrons/nm, because the slopes of the histograms in FIG. [3] differed; this was 
most likely due to background neutrons. 

Based on the £t observed between 2.0- and 6.0-nm wavelength for the data obtained 
with the supermirror guide, the most probable values of a and <|) were 0.448 ± 0.006 nm“^ 
and 2.0 x 10® neutrons/cm^/s, respectively. The value of T was fixed to 57.3 K, which was 
calculated via simulation of the moderator j^, 28|. Note that the flux around the 2.9-nm 
(136-m/s) components increased 30-fold compared with the ht results for the supermirror 
guide and B 4 G collimator. 


3.2. VCN transportation 


Six m = 3.5 wide-band multilayer monochromatic mirrors, fabricated at the Kyoto 
University Research Reactor Institute (KURRI), were installed in front of the unpolar¬ 
ized branch, as shown in FIG. 01 This was to ensure that only useful neutrons for UGN 
production were extracted, and also to decrease the background radiation. The reflectivi¬ 


ties of the monochromatic mirrors, as measured at SOFIA (BL16) at J-PARG 29|, l30|], are 


shown in FIG. O These values are plotted as a function of neutron momentum transfer 
Q = (dvr/A) sin0, where 6 is the glancing angle of an incident neutron. As shown in FIG. 01 
the direction of the VGN extraction is set to an angle of 18.6° below the horizontal. Hence, 
the mirrors are required to be arranged at an angle of (18.6° -|- 2 . 6°)/2 — 2.6° = 8.0°. Q 
is calculated as 0.80 nm“^ for a VGN with 2.9-nm wavelength making contact with the 
mirror at 0 = 10.6°, hence, the reflectivities are designed to be distributed around Q = 0.80. 
Accordingly, VGNs with approximately 2.9-nm wavelength are reflected downward, injected 
into the focus guide and then guided into the Doppler shifter, as will be described below. 

The mirror substrates in this apparatus are composed of 0.5-mm-thick Si with dimensions 
of 65 mm (width) x 65 mm (height) and are arranged with 11-mm spacing in order to cover 
the beam cross section within the shortest path. The mirror angle was hne-tuned in order to 
maximize the UGN yield, with the mirrors being tilted off the horizontal by 8.2°. Therefore, 
the beam makes contact with the monochromatic mirrors at a glancing angle of 10.7°. Note 
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FIG. 4: Arrangement of wide-band monochromatic mirrors. 
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FIG. 5: Reflectivity of wide-band monochromatic mirrors. 

that the mirrors are contained within a 5-cm-thick lead box covered by a 5-mm-thick B 4 C 
region, as shown in FIG. 121 to reduce the background noise from neutrons and gamma rays. 

As shown in FIG. El the VGNs reflected by the monochromatic mirrors are guided into the 
Doppler shifter by the Ni-coated guides, which consist of Ni-evaporated trapezoidal glasses. 
The guide cross section narrows towards the exit, so that the VGN flux is concentrated 
on the pseudo source point, while the divergence increases correspondingly. The guides are 
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composed of two parts: a long (515 mm) and short (230 mm) guide. The cross section of the 
VCN entrance is 60 mm (width) x 60 mm (height) and that of the exit is 40 mm (width) x 
26 mm (height). The exit size was optimized to maximize the UCN yield using a Monte 
Carlo simulation. The end of the short guide is positioned 182 mm from the pseudo source 
point. The TOF distance for the transit between the moderator surface and the pseudo 
source point is 17.2 m, therefore, a VCN with 136-m/s velocity arrives at the pseudo source 
point 127 ms after its generation. 

FIG. [7] shows the VCN flux at the exit window downstream of the Doppler shifter, 
as measured by the neutron beam monitor MNH10/4.2F with a square collimator of 
20 mm (width) x 30 mm (height) (see FIG. E]). The data sets were compensated based 
on the detector efficiency in order to convert the neutron ffux. 

The closed circles indicate the flux obtained under the present conditions, with the su¬ 
permirror, Ni-coated guides, and 6 wide-band monochromatic mirrors equipped, while the 
line without circles indicates the flux measured in the demonstration experiment 19(|, where 
only the B 4 C collimator and 4 narrow-band monochromatic mirrors were installed. The two 
plots are scaled to a beampower of 1 MW. We analyzed the last and succeeding time frames 
before the MR injections using the same approach described in section [TT] Then, the frame 
overlaps of the preceding pulse were subtracted, assuming that the preceding pulse had the 
same tail as the hnal pulse. 

The open circles in FIG. [7] indicate a VCN flux calculated with the VCN transport after 
the beam exit of the unpolarized branch using the UCN simulation. In this simulation, 
the initial x'- and y'-values [x' and y' are dehned in section 13.Ih were given uniformly, and 
the velocity divergences in the x'- and ^/'-directions were determined using the result of 
the pinhole measurement mentioned in section 13.11 The velocities in the branch direction 
were given according to the probability distribution in equation (| 6 ]), 0(A)/<h. The pulse 
shape and the delay by moderation of 0.10 ms listed in 28(| were taken into account in the 
UCN simulation. The monochromatic mirrors with the reflectivities shown in FIG. |5] and 
the focus guide were implemented. For comparison of the spectrum shape, the simulated 
count integrated between 2.5 and 6.0 nm was normalized to the measured count in the 
same region, where the overlaps by the preceding pulses were well removed. The simulated 
spectrum agreed well with the measurement. 

The 1 MW-equivalent simulated VCN fluxes 0^ with wavelength range of 2.8-3.0 nm. 
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FIG. 6: Doppler shifter and peripheral device layout. The inset is a macrograph of the wide-hand 
monochromatic mirror, focus guide, and UCN extraction guide arrangement. 
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FIG. 7: VCN flux reflected by monochromatic mirrors with and without upgrade, and simulated 
flux with upgrade, at Doppler shifter VCN exit window. Each data set is compensated based on 
each detector efficiency and scaled to a beampower of 1 MW. 

which is the target region for the Doppler shifter, at the position (i-iv) in FIG. [6] and 
measnred flnxes (pm at the position (i) and (iv) are snmmarized in TABLE. HI The meaning 
of each position nnmber is as follows: (i): the nnpolarized branch window, (ii): the exit of 
the focns gnide, (iii); the psendo source point, and (iv): the square collimator on the VCN 
exit window. The cross sections Sc used in the flux calculation at each position are also 
summarized for reference. The value of pm at the position (i) was estimated by a numerical 
integration of equation ([6]). tThe value of ps at the position (i) was normalized by pm at the 
same position. The decrease of Scps from the position (i) to (ii) was due to the the transport 
efficiency of the monochromatic mirrors and the focus guide, which were estimated to be 
64% and 84% for 2.9-nm VCN by the UCN simulation, respectively. The decrease of flux at 
the downstream of the position (ii) was due to beam diffusion in the free space of the Doppler 
shifter. In the estimation with the ratio of ps at the position (iii) to (iv), the measured flux 
at the position (iii) was found to be 8.7 x 10^ neutrons/cm^/s, which is 43-fold that of the 
demonstration experiment. The presumed flux, however, was half the simulated value at 
the same position. The loss is suspected to be due to the wide-band monochromatic mirrors 
and/or the focus guide. We are currently investigating the cause of the large flux decrease 
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observed in the experiment. 


TABLE I: Change of 1 MW-equivalent simulated VCN flux (pg and measured flux cprn with wave¬ 
length range of 2.8-3.0 nm in cross section Sc at the position (i-iv) in FIG. [6l ^The value of (ps at 
the position (i) was normalized by the measured one. 


Position 

Sc [cm^] 

(ps [cm ^s 

(pm [cm 1] 

(i) 

4.0 X 5.0 

t2.3 X 10® 

2.3 X 10® 

(h) 

4.0 X 2.6 

2.3 X 10® 

- 

(hi) 

3.0 X 3.0 

1.8 X 10® 

- 

(iv) 

2.0 X 3.0 

5.6 X 10^ 

2.8 X lO'^ 


3.3. Doppler shifter 


The structure and setup of the Doppler shifter are shown in FIG. [6l A blade is rotated 
like the hand of a clock by a servomotor in a cylindrical vacuum chamber. The vacuum 
around the motor shaft is maintained using magnetic fluid seals. A mirror holder is located 
at the tip of the blade and four neutron mirrors, which form the Doppler mirror, are stacked 
and mounted on the holder. The mirrors consist of d = 2.9-nm wide-band multilayers 


with 2.6-3.2-nm d-spacings 


16l-ll8| , and each mirror has a different reflectivity distribution. 


Each mirror has ion-beam-sputtered NiC/Ti multilayers on both surfaces of 0.3-mm-thick Si 
wafer [l6|, which was created at KURRI. The mirror size is 35 mm (width) x30 mm (height). 
Further, 2.5-mm-wide areas on both sides are covered by stoppers; hence, the effective 
reflection area of the Doppler mirror is 30 mm x 30 mm. 

FIG. |8] shows the reflectivities of the four Doppler mirrors as functions of wavelength 
measured at SOFIA, along with those of the stacked mirrors, which were estimated using a 
numerical calculation considering multiple reflections. The mirrors were set on the mirror 
holder from the bottom to surface in descending order of reflectivity. The reflectivity at 
Q = 2.16, which corresponds to 68.0 m/s, was 0.229 ± 0.004. 

The rotating radius of the mirror center is 325 mm and, in the experiment described here, 
the rotation frequency was set to 2, 000 rpm. The center of the mirror moves with a velocity 
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FIG. 8: Reflectivities of four m = 10 wide-band multilayer mirrors and numerically calculated 
reflectivities of the four stacked mirrors. 

of 68 m/s and the mirror interacts with incoming 25-Hz VCNs of 136-m/s velocity once in 
every three incident pnlses. Thus, the UCN production repetition frequency is 8.33 Hz. As 
a UCN of 3-m/s velocity travels a distance of 300 mm in 100 ms, this production frequency 
is suitable for the demonstration of UCN production with TOF measurements. 

The 12-MHz master clock of J-PARC and the 25-Hz continuous timing pulse are fed to 
a motor control unit and divided pulses are used to synchronize the mirror rotation with 
the VCN incidence. The phase offset of the rotation is controlled by the control unit in real 
time. The phase control jitter is maintained at less than ±3.5 mrad. 

The Doppler shifter has three UCN extraction ports, with one above and two beside 
the pseudo source point. As the UCN simulation showed that the decelerated neutrons are 
distributed on the parabola of equation ([5]), as shown in FIG. [Ib, the UCN detector was 
mounted on the upper port in this study, as shown in FIG. El The two side ports functioned 
as viewports for measurement of the rotating mirror timing using a pair of a laser and a 
photo-sensor, which indicated the timing when the mirror holder passed the pseudo source 
point. The accuracy of the timing detection was within ±0.1 ms. 

A UCN guide was inserted into the upper extraction port, which is connected from 
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26 mm above the center of the Doppler mirror to 30 mm in front of the UCN detector, in 
order to collect UCNs with a wide solid angle. The interior dimensions of the guide were 
30 mm (width) x 60 mm (depth) and the height was 143 mm. Ni-coated silicon wafers were 
glued on the guide inner surface. The TOF distance between the center of the Doppler 
mirror and the UCN detector window was 199 mm. The rotation plane was aligned along 
the center line of the unpolarized branch. The vacuum chamber can rotate around the motor 
axis from —5° to 30° to match the angle of incoming VCNs. The angle was set to 18.6° 
upward. During operation, the vacuum chamber was maintained at less than 3 Pa using a 
dry pump. 

3.4. UCN detector 

For UCN detection, we used a proportional gas counter, UCN DUNia-10, manufactured 
by A. V. Strelkov. The cylindrical gas chamber was hlled with ^He gas at 2.7 x 10^ Pa, and 
then with a gas mixture of 1% CH 4 and 99% Ar up to 1.1 x 10^ Pa. The detector is sensitive 
to slow neutrons only, as a result of the rarehed ^He gas. The chamber is 50 mm in length 
and has an entrance window of 90 mm in diameter placed at the bottom. The entrance 
window is covered with a pure A1 foil of 0 . 1 -mm thickness, and the foil is supported by 
metal beams on the outside of the chamber, for use in vacuum. In the chamber, a W anode 
wire is stretched in the radial direction, with a high applied voltage of 1.0 kV. Note that we 
replaced the built-in DUNia-10 preampliher with one that is insensitive to the microphonic 
noise caused by the mechanical vibrations of the motor rotation. Further, neutron output 
signals of +2 V in height and 5 /us in width are now input directly into a 16 ch input 
PHA-I-LIST system, 760-PRU06PIK, and recorded, along with the proton-delivery timing 
signals and the 1/3-divided 25-Hz pulses for the TOF analysis. 

The detection efficiency was evaluated by comparing the wavelength spectra measured for 
the cold neutrons by the UCN detector and those measured by a 1-MPa ^He proportional 
counter, RS-P4-0812-223, of 25.4-mm diameter with a 0.5-mm-thick stainless steel wall. 
Here, we label the attenuation by the A1 window and the stainless steel wall as ta and tb, 
respectively, and the detection efficiency of the UCN detector and the 1-MPa counter are 
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eA and e^, respectively. Then, the detection ratio of the two detectors, 'r(A), is given by 


r(A) = 


TAeA exp(-poA)(l - exp(-piA)) 


( 7 ) 


Tses exp(-p2A)(l - exp(-p3A)) ’ 

where po is the attenuation coefficient of the A1 window, pi is the reaction coefficient of 
the ^He gas of the UCN detector, p 2 is the attenuation coefficient of the stainless steel wall 
and ps is the reaction coefficient of the ^He gas of the 1-MPa counter. The values of po 


and p 2 were calculated to be 9.1 x 
A1 cross section at 298 K in Ref. 


0 and 6.0 x 10 ^ nni respectively, using the total 


3lj and the absorption cross section of Fe, instead of 


stainless steel, in Ref. |32|. Through htting in the 0.50-1.04-nni region, pi was evaluated as 
(1.001 ± 0.004) X 10“^ nm“^. The result was within 5% agreement with 9.5 x 10“^ nm“^, 
which was calculated from the ^He gas conditions of 2.7 x 10^ Pa and 50-mm thickness at 
300 K. Further, ps was estimated as 7.7 ± 0.6 nm“^. The detection efficiency decreased by 
20% when the incident point was 35-mm from the wire. This discrepancy was caused by 
insufficient collection of the electrons by the anode wire, as the pulse height was reduced by 
25% in the region far from the wire. 

As a result, the detection efficiency of the UCN detector without attenuation in the A1 
foil, e(A), can be expressed as 


e(A) = 1 — exp(—O.IOOA), 


( 8 ) 


with A in nanometers. 

4. MEASUREMENT OF DOPPLER-SHIFTED NEUTRONS 

In this section, we describe the results of UCN production experiments and the analysis 
using the UCN simulation. 

4.1. Doppler-shifted neutrons 

Typical TOF spectra for the VCN measured downstream of the operating Doppler shifter 
by the MNH10/4.2F neutron monitor are shown in the upper part of FIG. [9l Because the 
full time width of a VCN pulse separated by the wide-band monochromatic mirrors was 
approximately 100 ms, the spectra appear to be continuous because of the frame overlap. 
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As shown in the middle and lower parts of FIG. [9], the mirror rotates four times during 
the three VCN pulse incidences, so that four dips of 30-ms intervals appear every three 
40-ms pulses. These dips were used for a cross-check of the VCN reflection timing and an 
estimation of the reflected VCN velocity. The calculated error of the timing estimation was 
±0.2 ms. In the experiment, the phase offsets estimated using this method were consistent 
with those evaluated by the photo-sensor, as mentioned in section IT^ The timings obtained 
by these two methods agreed within the errors. 



FIG. 9: TOF spectra of VCNs passing through Doppler shifter (top). The four dips show the 
overlap of the Doppler mirror with the VCN beams. The middle and lower parts of this figure 
show the timing relation between the VCN pulse incidence and the VCN reflection by the Doppler 
mirror, schematically. 

The TOF spectra of the Doppler-shifted neutrons from the pseudo source point with 
phase offsets of 0.0, 2.0, 4.0, 5.0, 6.0, 8.0, and 10.0 ms are shown in FIG. dU) The timing 
origins of the horizontal axis for the TOF spectra were adjusted to each phase offset. The 
origin of the offset was determined at the timing when the photo-sensor pulse coincided with 
the 25/3-Hz timing signal. 

The TOF spectra with offsets of 0.0-6.0 ms have peaks at 15 ms. The peak height begins 
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FIG. 10: TOF spectra of Doppler-shifted neutrons with different phase offsets. Each time origin 
was set to the time when the Doppler mirror holder passed the pseudo source point. 

to decrease with increasing phase offset, and a peak at 45 ms arises. This is because the four 
dips in FIG. [9l which imply neutron reflection by the Doppler mirror, move to the right in 
response to the increased phase offset, as shown in the lower part of the hgure. Therefore, 
the neutron reflection timing gradually becomes mismatched with any given VCN pulse, but 
then begins to match with the next pulse. The timing relation between the dips and VCN 
pulses in FIG. |9] shows that the spectra with 0.0- and 10.0-ms offsets differ in this way only, 
but are otherwise intrinsically identical. Note that, in these measurements, the length of 
the UCN extraction guide was not 143 mm, but 121 mm. 

The neutron backgrounds were evaluated through measurement with the off-position of 
the Doppler mirror from the beamline and were then subtracted from each spectrum. The 
background count rate was 2.8 cps with a beampower equivalent to 1 MW, while that of the 
Doppler-shifted neutrons for the 5.0-ms offset was 4.1 x 10^ cps. Hence the background was 
negligibly small. The small peaks that appear at 60, 90, and 120 ms in the hgure represent 
the background neutrons scattered by the mirror holder; this can be concluded based on 
the fact that the small peaks appear immediately after the mirror crosses the VCN pulse. 
It is very difficult to block these background neutrons in the detector arrangement, because 
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they enter the UCN extraction guide with the UCN. Although the background neutrons can 
be evaluated precisely by utilizing plain silicon wafers instead of the Doppler mirror, the 
Doppler mirror can not be dismounted from the holder easily. 

The neutron yield between 31 and 35 ms, for which v was sufficiently less than that of 
the UCN and the spectrum did not overlap the peak from next pulse, was at a maximum 
for a phase offset of 5.0 ms. Therefore, the long-term experiments discussed below were 
performed with a 5.0-ms phase offset. 

The longitudinal wavelength spectrum of the Doppler-shifted neutrons with 5.0-ms phase 
offset is shown in FIG. [TTl along with that of the demonstration experiment. The wavelength 
was calculated based on the TOF distance and the duration between the offset and detection 
times, while considering the deceleration due to gravity. The flight path in the detector 
chamber was not considered in this figure. Each data set was compensated against the UCN 
detector efficiency given in equation ([8]), while the attenuation in the A1 window was ignored. 
The background values were subtracted from each plot. The small peak at 80 nm represents 
the abovementioned second peak arising from the reflection of the next VCN pulse. The 
peak that appears at 130 nm is due to the background neutrons scattered by the Doppler 
mirror holder. The 1 MW-equivalent total output was 6.0 x 10^ cps. Note that the phase 
offset of the demonstration experiment was estimated as 4.5 ms, from the VCN dips, and 
the 1 MW-equivalent total count rate was 13 cps. 

The part of the spectrum correspon ding to wavelengths longer than 123 nm, which is 
the cutoff wavelength of the A1 window 33f|, was regarded as a constant background. This 
background was most likely due to the neutrons remaining in the vacuum chamber, which 
would have been scattered by the rotating blade or the mirror holder and would have then 
entered the detector through the UCN extraction guide. These suspect neutrons surely 
mixed in the UCN region; thus, we assumed that a constant background with the same 
height as the 124 nm count increased the overall UCN spectrum. The constant background 
count between 58 and 123 nm was estimated as 7.8 cps. The total count of the second peak 
was estimated to be 2.1 cps; this estimation was conducted by approximating the peak as 
triangular with linear slope. Accordingly, the 1 MW-equivalent count rate for neutrons with 
longitudinal wavelengths between 58 and 123 nm was estimated to be 1.6 x 10^ cps. In 


the demonstration experiment 


19[ |. the neutron count calculated in the same manner was 


1.8 cps. Thus, the count rate increased 91-fold from demonstration experiment. 
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FIG. 11: Spectra of Doppler shifted neutrons. The closed circles show the result with the upgrades, 
while the crossed symbols show a demonstration experiment result (before the upgrade). 

4.2. Comparison with simulation 

The measured and simulated spectra of the Doppler-shifted neutrons at the UCN detec¬ 
tor position with 5.0-, 6.0-, and 7.0-ms phase offset are shown in parts (a-c) of FIG. [I2l 
respectively. The dashed lines indicate neutrons with wavelength longer than 58 nm (true 
UCN) calculated via the UCN simulation. The measured spectra were compensated based 
on the UCN detector efficiency given in equation ([8]) only. 

In the UCN simulation, the VCNs were seeded at the unpolarized branch window in the 
same manner as described in section [321 Only the VCNs with 120-155-m/s velocity were 
calculated, in order to reduce the computation time, under the condition that the hrst and 
second peaks were fully produced. The reflectivity of the Doppler mirror was taken from 
that of the stacked mirrors in FIG. [HI and the mirror size was set to 30 mm x 30 mm. 
Scattering of the VCNs by the mirror holder was ignored. The neutron path length in the 
detector gas chamber was calculated using equation (jH]), up to the detector length of 50 mm. 
The longitudinal wavelength was calculated by assuming the TOF distance was 199 mm, 
while the TOF was the sum of the flight times inside and outside the detector, as in FIG. 
EH The time origin of the TOF was set to the moment when the Doppler mirror passed 
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FIG. 12: Comparison of simulated and experimental results of longitudinal UCN wavelength spec¬ 
tra at phase offsets of (a-c) 5.0, 6.0, and 7.0 ms, respectively. The dotted, solid, and dashed lines 
indicate measured data, simulated spectra, and the simulated true UCN spectra, respectively. The 
UCN simulation results are normalized by the measured VCN amount between 2.5 and 6.0 nm. 

through the pseudo source point. The detection efficiency was calculated using the true 
neutron wavelength, which was computed from the magnitude of the velocity vector in the 
UCN simulation. One detection event was weighted by this efficiency and then corrected in 
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the same way as in the experimental data analysis. 

As can be seen in FIG. [12] (a), for the 5.0-ms offset, the simulated (solid line) and 
measured (dotted line) spectrum shapes are in good agreement for wavelengths less than 
94 nm. Disagreement occurs for wavelengths longer than this value, which is assumed to 
be caused by the metal beams supporting the foil. The measured count in the UCN region 
without the second peak, between 58 and 78 nm, was 6.5% less than the simulated one. By 
considering the velocity divergence of the UCN, the simulated true UCN count (dashed line) 
around 60 to 70 nm was significantly smaller than the measured longitudinal UCN count. 
This was because Doppler-shifted neutrons with wide divergence were collected by the UCN 
extraction guide and, hence, many neutrons around the longitudinal UCN limit were faster 
than 6.8 m/s. 

By normalizing the total count of the simulated spectrum by the measured result, the 
dashed line in FIG. [12] (a) is scaled to the true UCN spectrum in the experiment. In this case, 
the UCN count rate was estimated to be 92 cps. However, this value is slightly larger than 
the true UCN count in the experiment, because the measured spectrum had the constant 
background mentioned in section 14.11 and the spectrum around the cutoff wavelength of 
A1 was smaller than the simulated one. The constant background in the true UCN region 
and the disagreement around the cutoff wavelength were estimated to be 7.8 and 3.6 cps, 
respectively. After the removal of the counts from 92 cps, the true UCN count rate in the 
experiment was estimated to be 80 cps. This is half of the measured longitudinal UCN count. 
As a result, the UCN count rate increased 45-fold against the demonstration experiment. 
This is approximately the same as the VCN increase at the pseudo source point. 

The shape of the simulated spectra for the 6.0- and 7.0-ms offsets are also well matched 
to the experimental results, including the occurrence of the second peak, as shown in FIG. 
[12] (b) and (c). In the same manner as the 5.0-ms offset, the true UCN count rates for 
the 6.0- and 7.0-ms offsets were estimated to be 72 and 60 cps. The simulated hrst peak 
for the 7.0-ms offset case was 7.2% larger and the second peak was 24% smaller than the 
experiment. As the discrepancy was minimized for the 7.3-ms offset time setting in the 
UCN simulation, the TOF distance of the incident beam was 4 cm shorter or the reflectivity 
distribution of the Doppler mirror was 1% higher than our estimation. The systematic error 
of the phase offset measured by the photo-sensor could be also a source of the discrepancy. 
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5. UCN DENSITY ESTIMATION 


The maximum phase space density of the UCNs is limited to that of the incident beam, 
because a conservative force such as those acting in the Doppler shifter or the gravity does 
not cause a change in phase space density. Therefore, in this section, we estimate the phase 
space densities of the VCN beam and the produced UCNs, along with the number density 
of the UCNs at production in the same manner as Refs. 


10 
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The phase space density of the pulsed neutrons is calculated as follows: if a group of 
neutrons with various velocities is produced in a neutron source, the phase space distribution 
within a certain velocity range in the z-Vz plane is expressed simply as a trapezoid-shaped 
distribution (see FIG. ITT]) . Here, V 2 and Vi are upper and lower velocity limits, respectively, 
Z 2 and zi are the head and tail of the neutron pulse position, respectively, and tw is the time 
width of the neutron pulse. In the TOF measurement, the number of neutrons {N) between 
zi and Z 2 is ordinarily regarded as being those between vi and V 2 - 



FIG. 13: Schematic diagram of phase space distribution of pulsed beam in z-Vz plane. 

If the repetition frequency of the pulse is /, the duration of the beam per 1 s is ftw 
Therefore, the peak phase space density inside the trapezoid of the pulsed beam is 1/ft^, 
times greater than that of a continuous beam with the same average flux. The spatial 
volume of the beam passing through a unit area of the neutron detector per unit time is 
approximately described as the mean velocity v = (vi + V 2 )/ 2 , if (v 2 — vi) is sufficiently small. 
In this case, N represents the average flux of the pulsed beam. In addition, the divergences 
Avx of Vx and Avy of Vy are used to calculate the velocity phase space volume. Accordingly, 
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the peak phase space density, p(n), is expressed as 


p{v) = • V ■ Av^ ■ Avy ■ {v 2 - vi)). 


( 9 ) 


5.1. Phase space density of VCN beam 


The parameters used in the equation ([9]) and calculated peak phase space densities at 
the unpolarized branch window and the pseudo source point are summarized in TABLE. 
m The common parameters were estimated as follows. The pulse width (FWHM) was 
estimated for 136-rn/s VCN by htting the time spectra of the official moderator simulation 


for the BL05 port 


28| . The repetition frequency of the pulsed beam was as mentioned 


m 


section [31 The parameters at the unpolarized branch window were estimated as follows. The 
1 MW-equivalent flux between 134 m/s (2.95 nm) and 138 m/s (2.87 nm), was calculated 
from equation ([6]) without the constant Cb- The FWHMs of VCN divergences mentioned 
in section 13.11 were selected as the x'- and ?/'-divergence. The divergences were converted 
into the velocity widths of 136 m/s-VCN, namely Avx' and Avy'. The ratio of VCNs in 
the Avx'- and An^'-region to total VCNs in the velocity phase space was named “selected 
neutron yield”. The flux multiplied by the selected neutron yield was substituted for N in 
equation ([3]). The parameters at the pseudo source point were estimated as follows. The 
1 MW-equivalent flux that passes the cross section of 30 mm x 30 mm was estimated from 
the UCN simulation and the data set at the square collimator on the VCN exit window 
as mentioned in section 13.21 The x- and ^/-divergence of VCN were determined by the 
boundaries of dense region in histograms of calculated beam divergence. The divergences 
were converted into Av^ and Avy. The meaning of “selected neutron yield” and its usage 
are the same as above mentioned. 

As shown in TABLE. ini the phase space density at the pseudo source point is 15% of 
that at the unpolarized branch. However, the value is more than 3 times larger than the 
UCN phase space density of 8.4 x 10“^ neutrons/cm^/(m/s)^ at ILL/PF2 j^]. If the flux 
at the unpolarized branch window is ideally transported to the pseudo source point, the 
efficiency is calculated to be 76% from TABLE. HI Hence, the flux at the pseudo source point 
was expected to be 7.4 x 10'^ neutrons/cm^/s. However, the measured flux, and accordingly 
measured phase space density, was half the expected value. 

Taking the reflectivity of the Doppler mirror as being 0.229, the phase space density of 
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TABLE II: Parameters used in the calculation of peak phase space density and the result of 
the calculation at the unpolarized branch window and the pseudo source point, x'-, y'-, x- and 
y-divergence were converted into the velocity widths of Avx', Avy/, Avx and Avy, respectively. 
“Selected neutron yield” means the ratio of VCNs within the Avx'- (or Avx-) and Avy/- (or Avy-) 
region to total VCNs in the velocity phase space. 


Common parameters 

Pulse width (FWHM) 

4.9 X 10^ ys 

Repetition frequency 

25 Hz 

Parameters at the unpolarized branch window 

Neutron flux in 134-138 m/s 

9.6 X 10^ cm“^s“^ 

x'-divergence (FWHM) 

50 mrad 

y'-divergence (FWHM) 

54 mrad 

Avx' 

6.8 m/s 

AVy' 

7.3 m/s 

Selected neutron yield 

61% 

Phase space density 

1.8 cm“^(m/s)“^ 

Parameters at the pseudo source point 

Neutron flux in 134-138 m/s 

3.5 X 10^ cm“^s“^ 

x-divergence 

0 ± 24 mrad 

y-divergence 

0 ± 70 mrad 

Avx 

6.5 m/s 

AVy 

19 m/s 

Selected neutron yield 

63% 

Phase space density 

0.26 cm“^(m/s)“^ 


the reflected VCN is expected to be 6.0 x 10“^ neutrons/cm^/(m/s)^. The value is nearly 
the same as the phase space density of UCN at ILL/PF2. If the UCNs hll the entire 
spherical velocity volume with 6.8-m/s radius, the UCN number density is calculated as 
being 80 UCN/cnU. 
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5.2. UCN phase space and number densities 


The phase space density of the detected UCN was estimated as follows. The 1 MW- 
equivalent measured UCN count was estimated to be 9.6 UCN/pulse from the conclusion 
of section 021 The spatial volume where the UCNs were produced was 3.0 cm x 3.0 cm x 
6.0 cm according to the UCN simulation. The velocity phase space in which the UCNs 
were produced has a partial spherical shell shape with radius from 3.7-7.1 m/s, determined 
by considering the A1 cutoff and gravity. The volume of the phase space was numerically 
calculated to be 163 (m/s)^. Therefore, the phase space density of the UCN at production 
is 1.1 X 10“^ neutrons/cm^/(m/s)^, which is 1.8% of the phase space density of the VCN 
reflected by the Doppler mirror. This corresponds to 1.4 UCN/cnU, assuming a spherical 
velocity phase space volume of 6.8-m/s radius. The density is 12 times larger than that 
produced at the Argonne National Laboratory [h| and is as high as that in the nEDA 
experiment at ILL/PF2 which set the current upper limit of 2.9 x 10“^® e • cm (90% C.L. 

3 - 


The source of the reduction in the phase space density can be explained by an overesti¬ 
mation of the phase space volume because of empty components contained within it. During 
the UCN production, the Doppler mirror overlaps the injected VCN beam for only 0.9 ms, 
as mentioned in section [2l so the range of the reflected VCN velocity in the z-direction is 
1 m/s. The UCN range is identical to that of the VCN expressed as a black parallelogram 
in FIG. [131 if we assume a Doppler mirror moving parallel to the VCN for simplicity. 

In the UCN simulation, the FWHM of the distribution of the reflected VCNs was 
0.58 m/s. Therefore, the velocity phase space occupied by the UCN has a disk shape with 
0.58-m/s thickness and 6.8-m/s radius. The volume of the disk is 84.9 (m/s)^, which is 6.4% 
that of a sphere with radius 6.8 m/s. By multiplying the VCN phase space density obtained 
in section l5T] and the velocity volume of the disk together, the UCN density is calculated as 
5.1 UCN/cm^. The spatial volume was also overestimated because the UCN density in the 
volume was not uniform. In the UCN simulation, the number density of UCNs in the 1-cm^ 
cube at the psuedo source point just after production was 3.9 UCN/cm^, while the average 
density in 3.0 cm x 3.0 cm x 6.0 cm-volume was 1.2 UCN/cm^. Thus, the overestimations 
in the spatial and velocity volumes roughly explain the reduction of the phase space density 
of the UCNs. 
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As the UCNs produced by the Doppler shifter have a position and velocity correlation, 
they are concentrated in a part of the phase space volume. As described in equation (E]), 
the Doppler-shifted neutrons are distributed on a line of {2vm — v) in the y-v^ plane, as 
shown in FIG. [TU because of the difference of the rotating radius. Therefore, if we take the 
phase space volume as a product of the spatial and velocity volume with which the UCNs 
are produced, the effective phase space density becomes small. If we use a parallel moving 


Doppler mirror, such as an excentric double motor 


34l |. or modify the UCN velocity with 


a magnetic accelerator system like the UCN rebuncher, we can obtain a UCN beam with 
smaller divergence. 
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FIG. 14: UCN distribution in y-Vz plane as calculated by the UCN simulation. The dashed line 
shows the dominant term of equation ([5|) in the region, that is {2vm — v) where v = 17.2/0.125 ~ 
138 m/s and Vm = (2007r/3)y. 

6. CONCLUSION AND FUTURE OUTLOOK 

We have developed the proposed Doppler shifter for the hrst UCN source at J-PARC and 
signihcantly improved the UCN output by replacing or adding neutron optical components 
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for practical use. The Doppler-shifted neutrons have an 8.33-Hz sharply pulsed time struc¬ 
ture of 4.4-ms width. The 1 MW-equivalent neutron count rate for neutrons of longitudinal 
wavelength between 58 and 123 nm is 1.6 x 10^ cps. The true UCN count rate is estimated 
to be 80 cps, according to a Monte Carlo simulation. The UCN number density is estimated 
to be at least 1.4 UCN/cm^ at production. The Doppler shifter works almost as designed 
and produces sufficient UCN output for practical use. 


Through application of the Dopp 
for use in our search for the nEDM 


er shifter, we have already begun R & D of apparatus 


20[ |. The result of the UCN transport experiment with 


the 5.6 m-long guide is shown in FIG.dSl as an example. UCNs were guided to the horizontal 
direction with an L-shape mirror guide connected to the 5.6 m-long guide. UCN count was 
reduced to 1/3 in this section. The mean reflectivity of the 5.6 m-long guide was estimated 
to be 97%. The 1 MW-equivalent total count rate increased to 1.9 cps by the improvement 
of the Doppler shifter and the UCN region was 0.39 cps, although the UCN pulses were 
thinned to 1/16, 0.52 Hz, with the shutter. As the result of the improvement, the Doppler 
shifter allows us to perform the UCN rebuncher experiment at BL05 in a practical time 
span. At present, we are planning the performance test of the improved UCN rebuncher 
with an experimental setup incorporating the UCN guides, along with the development of 
a UCN reflectometer and storage experiments. 


7. ACKNOWLEDGEMENTS 

This research was partially supported by Ministry of Education, Science, Sports, and Cul¬ 
ture Grants-in-Aid for Scientihc Research (A) 23244047 and (A) 26247035, and by the JSPS, 
Grant No. 23360428. The Doppler shifter experiment was approved by the Neutron Scat¬ 
tering Program Advisory Gommittee of IMSS, KEK (Proposal No. 2009S03 and 2014S03). 
The neutron reflectivity experiment was approved by the Neutron Science Proposal Review 
Gommittee of J-PARG/MLF (Proposal No. 2009S08 and 2014B0263) and supported by the 
Inter-University Research Program on Neutron Scattering of IMSS, KEK. 


30 






TOF (sec) 



FIG. 15: TOF spectrum of UCN transport experiment with Ni-coated guides of 5.6 m in length. 
The UCN pulses were thinned to 0.52 Hz. 
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